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ABSTRACT 

We compute the hard X-ray spectra from a hot plasma pervaded by small cold dense 
clouds. The main cooling mechanism of the plasma is Compton cooling by the soft 
thermal emission from the clouds. We compute numerically the equilibrium tempera- 
ture of the plasma together with the escaping spectrum. The spectrum depends mainly 
on the amount of cold clouds filling the hot phase. The clouds covering factor is con- 
strained to be low in order to produce spectra similar to those observed in Seyfert 
galaxies and X-ray binaries, implying that an external reflector is required in order 
to reproduce the full range of observed reflection amplitudes. We also derive analyt- 
ical estimates for the X-ray spectral slope and reflection amplitude using an escape 
probability formalism. 
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1 INTRODUCTION 

The physical conditions in the inner parts of the accretion 
flow surrounding a black hole are likely to be very chaotic. 
A situation that has been often considered in the literature 
is that of the so-called 'cauldron', where a soup formed by a 
hot plasma contains small grains constituted by small dense 
clouds of much colder matter (e.g. Guilbert & Rees 1988). 
Several works were devoted to explain how such a configura- 
tion could be physically realized and compute the spectrum 
emitted by the clouds for different cloud optical depths (see 
e.g. Rees 1987; Ferland & Rees 1988; Rees, Netzer & Fer- 
land 1989; Celotti, Fabian & Rees 1992; Kuncic, Blackman 
& Rees 1996; Collin-Souffrin et al. 1996; Kuncic, Celotti & 
Rees 1997; Krolik 1998). In particular the main observable 
effect of optically thick cold clouds is the reprocessing into 
soft UV photons of the hard X-ray radiation produced in the 
hot phase by the Comptonisation process. This reprocessed 
emission is likely to contribute, at least partly, to the big 
blue bump observed in AGN. In addition, the presence of 
such clumps inside the hot plasma may also contribute to 
the formation of a reflection component (see e.g. Nandra & 
George 1994), responsible for the bump in the hard X-ray 
domain, commonly observed in Seyfert galaxies and galactic 
black hole candidates. 

Most of the previous works on cloud models focused on 
the physics and radiative processes in the clouds themselves, 



without taking into account the possible effects of the clouds 
on the characteristics of the hot phase. Indeed the soft ra- 
diation re-emitted by the cold material can constitute the 
main radiation field responsible for the Compton cooling of 
the hot gas. This in turn affects the temperature of the hot 
plasma, and thus the emitted X-ray Comptonised spectrum. 

This feedback loop is identical to that found in accre- 
tion disc corona models (Haardt & Maraschi 1993). Simi- 
larly, the conditions at radiative equilibrium depend mainly 
on the cold matter distribution relative to the hot plasma. In 
this context, Malzac (2001) (hereafter M01) studied a geom- 
etry where the clouds are external to the hot Comptonising 
plasma and spherically distributed around it. The aim of the 
present work is to study the effects of the presence of cold 
optically thick clouds distributed inside the hot phase. As 
demonstrated below, the cooling by the cold clouds is then 
more efficient than in the case of an external reprocessor. 

Under assumptions detailed in section ^, we use a nu- 
merical approach based on non-linear Monte-Carlo simula- 
tions, described in section ^, to compute the emitted spectra. 
These numerical results are found to be in agreement with 
analytical formulae, derived in section ^, giving the slope 
of the primary X-ray Comptonised spectrum V and the re- 
flection amplitude R. The predictions of the model are then 
discussed and compared with the data in section [B|. 
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Figure 1. The hot plasma emitting the hard X- and soft 7-ray 
emission is pervaded by small lumps of cold dense matter. The 
geometry of the hot plasma could be close to spherical (left) or 
slab (right). The cold matter is assumed to be homogeneously 
distributed inside the plasma volume. The ambient high energy 
radiation which is intercepted by the cold clouds is partly repro- 
cessed as low energy (UV, EUV) radiation and partly reflected in 
the X- and 7-ray energy domains. The system is assumed to be 
in radiative equilibrium. 



2 MODEL ASSUMPTIONS AND 
PARAMETERS 

The physical characteristics of the clouds are of prime impor- 
tance for the emitted spectrum. If they are optically (Thom- 
son) thin, and with a large covering fraction, partial trans- 
mission of the X-ray radiation across a cloud is likely to 
produce strong absorption features in the X-ray spectrum 
(e.g. Kuncic, Celotti & Rees 1997). On the other hand, if 
the clouds are optically thick and highly absorbing they do 
not produce, independent of their covering fraction, strong 
apparent absorption features, as indeed observed (Nandra 
& George 1994). 

For this reason in this analysis we will consider the case 
of optically thick clouds with zero transmission, which corre- 
sponds to typical column densities <; 10 25 cm -2 . We assume 
that the individual cold clouds are much smaller than the 
characteristic size of the region occupied by the hot plasma: 
if the emitting region has a typical dimension H, the typ- 
ical cloud size is of order eH with e & 10~ 2 ^. This con- 
strains the cloud hydrogen density to be large, larger than 
(oieH) -1 . This limit writes uh <; 10 12 cm~ 3 in the case 
of Seyfert galaxies (AGN, with H ~ 10 15 #15 cm), and 
uh il 10 20 cm -3 in the case of galactic black holes (GBH, 
where H ~ 10 7 Hr cm). Such large densities in turn con- 
strain the ionization parameter to be relatively low (£<10 3 
for a 10 45 erg s _1 luminosity AGN, and £<10 2 for a 10 36 
erg s~ GBH) and thus the possible effects of ionization are 
likely to be weak (Zycki et al. 1994). In our calculations we 
will then assume that the clouds are neutral. 

We will consider a system in radiative equilibrium. The 
effective temperature of the clouds fcTbb can be simply esti- 
mated. Assuming that the flux is uniform inside a spherical 
hot phase of bolometric luminosity L, a cloud surface el- 
ement will receive a flux F = L/(4nH 2 ). If the clouds arc 
dense enough this flux is fully absorbed and re-radiated as a 
quasi-blackbody (e.g. Ferland & Rees 1988) and the Stefan- 
Boltzmann's law gives well known characteristic energies: 



* For simplicity we will assume spherical clumps, although fila- 
mentary or sheet-like geometries might be also appropriate. 
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where L45 and L36 are the bolometric luminosity expressed 
respectively in 10 45 and 10 36 erg s _1 units. We will assume 
that the radiation reprocessed by the clouds is the only 
source of soft seed photons. 

As we are focusing on the spectral signatures of this 
two phase medium and as they do not constrain - as al- 
ready mentioned - the physical status of the clouds (as long 
as they are optically thick and dense enough to maintain 
a low temperature), the clumps might exist with a large 
range of conditions. Nonetheless the range of parameters 
considered (in dimension and density) are consistent with 
assuming that the gas reaches radiative equilibrium (i.e. the 
cooling timescales are shorter than the dynamical one) and 
pressure equilibrium (with a size eH small enough to reach 
it in a dynamical timescale). 

In general if neutral clouds are confined by a hot ionized 
environment the thermal pressure balance can be written in 
terms of the cloud Thomson optical depth r = riH<JTeH: 



t ~ 2er T T e /Tb' 



(2) 



For typical observationally inferred values of Tf = 1 
and fcT e =100 keV, and using the expressions (^) for the 
temperature of the cold clouds, equation ^ reads: 



64e_ 



H ' 



1/2 



-r/4 



12e- 2 



1/2 



- 1/4 



(3) 



where e = 10~ e-2. We note however that this simple esti- 
mate does not consider the possibility of a two temperature 
plasma and the consequent effects of hot ions on the thermal 
balance (see Zdziarski 1998; Spruit & Haardt 2000; Deufel 
& Spruit 2000). 

Thus the thick clouds considered here can be confined 
by the hot phase, but in any case thicker clouds may ex- 
ist and be e.g. magnetically confined (for more details see 
e.g. Kuncic, Blackman & Rees 1996). Note that the clumps 
do not need to survive for longer that their cooling time 
provided they can be continuously replaced. 

Let us now consider the spectrum produced in these 
environment. The three main parameters controlling the re- 
sulting X-/7~ray spectral shape are the following: 

(i) The amount of cold clouds pervading the hot plasma. 
Since we consider small scale clumps we can quantify the 
amount of cold matter by using an effective "cloud" optical 
depth tb, such that a photon crossing the hot plasma has a 
probability 1 — exp(— tb) of intercepting a cloud, tb may be 
thus defined using a cross section formalism, considering an 
individual cloud as a particle: 



N 

r B = y (A)H 



(4) 



where V is the plasma volume, TV is the total number of cold 
clouds, (A) is the average geometric cross section of each 
cloud. For spherical clouds pervading a spherical plasma, tb 
can be written explicitly in term of the number and size of 
the clouds: 



TB = 



3/Ve 



(5) 



(ii) The Thomson optical depth tt of the hot plasma. 
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Figure 3. The volume averaged electron temperature versus tb. 
Filled and open symbols correspond to fcTbb = 150 eV (GBHs) 
and fcTbb = 5 eV (AGN), respectively. Squares and triangles refer 
to the cases tt = 2 and tt = 0.5. 



(iii) The characteristic energy of the soft photons emit- 
ted by the clouds. We assume a blackbody spectrum with a 
fixed temperature fcTbb. As the radiative equilibrium is not 
extremely sensitive to the value of fcTbb, we will fix fcTbb=5 
eV and fcTbb = 150 eV as typical values for Seyfert galaxies 
and GBH sources, respectively, as from equation (Q). 

As said before, the amount of cold matter and its dis- 
tribution with respect to the hot plasma are of prime im- 
portance. On the other hand, the actual geometry of the 
emitting region has little effects on the emitted spectrum. 
We studied both the spherical and infinite slab geometries 
illustrated in Fig. |l|. As the results were qualitatively sim- 
ilar for both geometries, in the following we consider the 
spherical case only. 



3 NUMERICAL APPROACH 

In order to estimate the resulting physical parameters and 
spectra we use the Monte-Carlo code of Malzac & Jour- 
dain (2000) which is based on the Non-Linear Monte-Carlo 
method (NLMC) proposed by Stern (Stern 1985; Stern et 
al. 1995). We consider a spherical geometry. The sphere is 
divided in 5 shells with equal thicknesses and is assumed to 
have a uniform density. Both tb and tt are defined along 
the radius H of the sphere. 

In the hot phase, the simulations are fully non-linear: 
both photons and electrons are followed using the Large Par- 
ticles (LP) formalism described in Stern et al. (1995). This 
enables the equilibrium temperature to be computed in each 
zone according to the local energy balance heating = Comp- 
ton cooling. 

The interaction between a photon and a cloud is dealt 
as a regular Monte-Carlo interaction between two particles. 
However, as the cold clouds structure is assumed to be un- 
affected by the radiation field, the clouds themselves are not 
followed. As the clouds are small and homogeneously dis- 
tributed, tb is simply the reaction rate for 1 photon, with 



time expressed in units of the light crossing time H/c of the 
medium. When a LP photon interacts with a cloud, it is 
assumed to enter a semi-infinite slab medium with standard 
abundances and neutral matter (Morrison & Mc Cammon 
1983). The linear Monte-Carlo code of Malzac et al. (1998), 
is then used to track the path and the energy changes of 
the LP photon in this medium, until it is either absorbed or 
escapes from the cold matter slab. As long as the clouds are 
optically thick, the reflection spectrum on a cloud does not 
differ much from slab reflection (Nandra & George 1994). A 
more detailed modeling of the reflection component would 
rely on a specific cloud geometry which is essentially un- 
known. 

The energy deposited in the cloud is re-injected in the 
hot phase in the form of thermal soft photons LP. Both the 
reflected and reprocessed LP have a direction of propagation 
drawn from an isotropic distribution. 



4 ANALYTICAL ESTIMATES OF R AND V 

There are two parameters generally used to phenomeno- 
logically describe the X-ray spectral properties of accret- 
ing black hole sources, namely the amplitude R of the re- 
flection component with respect to the primary continuum 
and the photon spectral index F of the power-law shaped 
primary component (see Malzac, Beloborodov & Poutanen 
2001, MBP hereafter). 

These parameters were determined from the numeri- 
cal simulations in a way similar to that described in MBP. 
The reflection was estimated by computing the flux ratio be- 
tween our simulated reflection, and the angle-averaged slab- 
reflection produced by a source having the same primary 
spectrum as the simulated one. 

However, using the simple escape probability formal- 
ism described in the following, it is also possible to indepen- 
dently derive analytical estimates for R and T, as follows. 

A power Lh is deposited in the plasma as heating and 
escapes the system as three different kinds of radiation: 



Lh — Luv + Lx + Lr, 



(6) 



respectively the UV/soft X-ray reprocessed luminosity, the 
X-ray (Comptonized) luminosity, and the reflected lumi- 
nosity. Let Uh be the total radiative energy inside the hot 
plasma. With similar notations as in equation (n) we have: 



u h = U UV + U X + U R . 



(7) 



The escaping luminosity is related to the internal en- 
ergy through the escape probability (Lightman & Zdziarski 
1987): 



P h = (L h /U h )(H/c) 



(8) 



which depends on the geometry, sources distribution, optical 
depth and energy. We consistently indicate as P x , P s ,Pr 
the escape probability for Comptonised, soft and reflected 
radiation and will use the analytical approximation to P 
given in the Appendix (equation (po|)). 

Let us now estimate the three different luminosities. 

The soft luminosity in the hot phase is produced by the 
clouds through absorption and reprocessing of the Comp- 
tonized and reflected radiation and disappears through es- 
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Figure 2. Examples of spectra produced by the lumpy model as a function of the indicated values of tb- The dashed lines show the 
resulting spectra without the reflection component. The left panel corresponds to the case of AGN (fcTbb = 5 eV). The right panel 
corresponds to GBHs (kT^b = 150 eV). The spectra (in arbitrary flux units) are angle-averaged and normalized to the same total 
luminosity. The plasma Thomson optical depth is tt = 2. 




Figure 4. Angle averaged reflection (left-hand panel) and spectral index (right-hand panel) versus tb - Filled and open symbols correspond 
to fcT bb = 150 eV (GBHs) and fcT bb = 5 eV (AGN), respectively. Squares and triangles show the cases tx = 2 and -pp = 0.5. In both 
panels, the curves represent the analytical approximation to the numerical results (see text). The solid lines stand for fcT bb = 5 eV 
(AGN), the dashed lines for kT hh = 150 eV. 



cape or Comptonisation. Thus, the radiative equilibrium 
balance for U uv reads: 

(P s + tt)[/ uv = Tb(1 - a)(7 x + tb(1 - a R )U R , (9) 

where a and an are the clouds energy and angle integrated 
albedo for a Comptonised and a reflection spectrum, respec- 
tively. For neutral matter the typical values are a ~ 0.1 
(Magdziarz & Zdziarski 1995; MBP) and a R ~ 0.4 (Malzac 
2001). 

Similarly, the reflected radiation is formed through re- 
flection of the Comptonised and reflected radiation on the 



clouds, and disappears via escape, Comptonisation in the 
hot plasma and absorption by the clouds: 

(Pr +t t + tb(1 - a R ))U R = r B aU x . (10) 

The fraction of the internal energy in the form of Comp- 
tonised radiation has as main source the power dissipated in 
the hot plasma but also the Comptonised soft and reflected 
radiation; the sink term is due to escape, absorption and 
reflection by the clouds: 

(A + r B )C/ x = L b H/c + r T U uv + t t Ur. (11) 
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Excluding the reflection component, the total luminos- 
ity sinking out of the hot phase either through escape or 
interaction with the clouds is: 



L = (P s + tb)U s + (P s + r B )U uv , 



(12) 



while the soft luminosity entering the hot phase (coming 
from the clouds) is given by equation (Eh: 



(Px + T T + T B )U U 



(13) 



Therefore combining equations ([[2]) and ( |l3| ) with equa- 
tions (^) and (^) provides an estimate for the amplification 
factor A = L/L s : 



A 



(P s + t b )t b + (P x + r B )(P s + r T ) 



Tb(P S + TT + Tb) 

where the factor g is given by: 
Pr + tt + tb(1 — an) 



9 = 



(14) 



(15) 



(Pr + t t )(1 - a) + tb(1 - a R ) 

Note that at first order g ~ 1 and P x — Ps and thus: 

A~l + P x /r B . (16) 

Finally the amplification factor is clearly directly re- 
lated to the photon index V, which can therefore be ex- 
pressed - using the phenomenological formula given by Be- 
loborodov (1999) - as 



T = 2.33(4 - 1) 



-1/6 



(17) 



where the parameter 8 = 1 / 10 for AGN and 1/6 for GBHs. 

In addition equation (lid) provides an estimate for the 
second parameter, the reflection coefficient: 



R 



Pr 



TB 



aL x 



Px Pr + tt + tb(1 - an.) 



(18) 



This formula corresponds to an angle averaged reference 
slab-reflection spectrum. Actually, the reflection coefficient 
also depends on the assumed inclination angle for the infinite 
slab model (see e.g. pexrav model in XSPEC, Magdziarz & 
Zdziarski 1995). Equation (|l|) could be corrected for this, 
for example by dividing it by the angular factor given by 
equation 2 of Ghisellini, Haardt & Matt (1994). However, 
for the relatively low inclination angles usually assumed in 
spectral fits, the correction is small and, for simplicity, we 
will neglect it. 



5 RESULTS 

Figure ^ shows some of the simulated spectra for different 
values of tb. The radiative cooling from the clouds appears 
to be extremely efficient and starts being effective for values 
of tb as low as a few x 1CP 2 . At fixed tt and Tbb, increasing 
tb increases the cooling of the plasma (because of the en- 
hancement of the reprocessing and the subsequent increase 
of the soft radiation field) and therefore decreases its tem- 
perature (see below). This affects the shape of the emitted 
spectrum which becomes softer and cuts off at lower en- 
ergies. Note that for the same set of parameters the spec- 
tra of GBHs (right panel of the figure) are harder than the 
AGNs ones (left panel) due to the larger energy domain 
where the luminosity can be released in the latter objects 



(see e.g. MBP). The spectra shown are rather similar to 
what is observed in Seyfert galaxies and Galactic black hole 
candidates in their hard state (see e.g. Zdziarski et al. 1997, 
Poutanen 1998). 

The dependence of the hot plasma temperature on tb 
is illustrated more precisely by Fig. [| As just said the tem- 
perature is generally lower in AGNs (for the same tb) and 
it is higher for lower Thomson optically depths. While these 
are general features (see e.g. MBP) of the self-regulated 
Compton emission generally invoked by the two-phase mod- 
els (Haardt & Maraschi 1993), note that in our case the 
product T c tt is not expected to be constant at fixed tb 
since here the Compton parameter depends also slightly on 
tt (see Section ^) . 

The observed range of temperatures - say 50-200 keV 
- inferred from observations of the high energy cut off in 
Seyfert galaxies and GBHs can be reproduced for tb in the 
range 10~ 2 -1, depending also widely on tt. In the optically 
thick case larger values of tb may induce very low temper- 
atures. 

A second observable parameter is the index of the (pri- 
mary) X-ray spectrum. Fig. ^ shows the evolution of the 
2-10 keV spectral slope for increasing tb: as expected the 
spectrum softens very quickly for larger tb . One can see that 
the range of spectral slopes observed in Seyfert galaxies and 
GBH in the hard state (approximately the range 1.4-2.2) 
can be achieved for tb roughly in the interval 10 _3 -1 (larger 
tb clearly lead to a too soft spectrum). Note that although 
such interval is rather large, the dependence on tb is quite 
steep especially in the GBH case. 

Finally let us examine the results on the intensity of the 
reflection component, shown in Fig. ^ (left-hand panel). This 
immediately reveals that the inferred values of R are rather 
low, even when tb approaches unity, and in general signif- 
icantly lower than observed. According to equation (|l8|), 
R values of order of unity can be achieved in the limit of 
large tb. However, this corresponds to very large (and un- 
observed) r, as discussed in the next section. 

As shown in the same Fig. ^, we also stress here that 
there is good agreement between the numerical results and 
the analytical estimates of R and T derived in section ^j. 
Concerning T,we note however a significant discrepancy for 
the GBH optically thin case, for which the simulations give 
a much harder spectrum than the analytical estimates. In- 
deed, in this case the first Compton scattering order falls in 
the 2-10 keV range. This forms a bump making the spec- 
trum locally harder in the energy range used to estimate F, 
and with a break at higher energy. This (unobserved) radia- 
tive transfer effect is due to the low optical depth (tt=0.5) 
considered. A part from this the main discrepancies appear 
in the optically thick case and at large tb (corresponding to 
r outside the observed range) arising because in this inter- 
val of parameters associated with low temperatures, equa- 
tion (p"7|), which is formally valid for a spherical plasma with 
temperature in the range 50-100 keV, does not provide a 
good approximation to the r vs A relation anymore. 



6 CONCLUSIONS 

From the computation of the spectra expected in a scenario 
where cold reflecting and reprocessing gas is embedded in 
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Figure 5. The circles show data in the _R~r plane corresponding 
to the sample of Seyfert galaxies (open circles) of Zdziarski et 
al. (1999) and the sample of galactic black holes (filled circles) 
presented in Gilfanov et al. (2000). The lines refer to the analyt- 
ical model predictions: along these lines tb is varied at constant 
Tt- Dashed and solid lines correspond to fcT bb =150 eV (GBHs) 
and fcT b b=5 eV (AGN) respectively. For both types of sources, 
the model predictions are plotted for tt=0.5 (upper curves), and 
tt=2 (lower curves), tb covers the range 0—2 along the curves. 



a coronal hot plasma we have inferred constraints on the 
distribution of such cold gas. In fact, having assumed that 
the gas can be described as small optically thick clouds, the 
main parameter regulating the amount of reprocessing (in 
turn responsible for the Compton cooling) and reflection is 
the cloud optical depth tb, which quantifies the covering 
factor of the cold component. 

This has to be relatively low. In fact, the observed range 
of spectral slopes constrains tb < 1. This limit however for- 
bids to produce simultaneously a reflection amplitude larger 
than ~ 0.3. The situation is even more critical since R ~ 0.3 
is predicted only for the softest sources which observation- 
ally appear to show the largest reflection amplitude (R^ 1), 
as shown in Fig. [s| 

Thus, in order to reproduce the spectral characteris- 
tics of numerous objects (both Seyfert galaxies and GBH 
sources) with significant reflection the model requires an ad- 
ditional reflecting medium which does not contribute to the 
soft photon field, such as e.g. a cold outer accretion disc, 
reflection on a torus in the case of Seyfert galaxies or on the 
companion star in the case of X-ray binaries. 

The reported correlation between R and V (shown in 
Fig. ^) is generally interpreted as being due to plasma cool- 
ing on soft photons emitted by the same medium that gives 
rise to the reflection component (see Zdziarski et al. 1999; 
Gilfanov et al. 2000). If the R-T correlation were indicative of 
a common medium for soft photon creation and reflection, 
such a medium cannot be provided by cold matter mixed 
together with the hot phase, and another reflector should 
contribute to, and even dominate the radiative cooling of 
the hot plasma. 

A related problem is the observed correlation between 



the degree of relativistic broadening of the Fe line with F 
(Gilfanov et al. 2000 for black hole binaries, Lubihski & 
Zdziarski 2001 for Seyferts, see however Yaqoob et al. 2002). 
In the context of the present model this correlation would 
suggest that the overall cloud distribution is not homoge- 
neous. The changes in the amount of cold material pervad- 
ing the hot plasma should be stronger in the innermost parts 
of the accretion flow where the relativistic broadening is im- 
portant. However, since the reflection amplitude produced 
by the cold clump is weaker than observed, this requires that 
the bulk of the line would be emitted by the cold clouds while 
the reflection would come mainly from an external reflector. 

Thus, the origin of both soft photons and reflection are 
unlikely to be matter mixed inside the hot Comptonising 
plasma. In the framework of the cloud models this suggests 
that the cold clouds should be essentially external to the hot 
Comptonising plasma. Then the radiative cooling is much 
less efficient and the reflection component may be larger. 
Variations in the covering factor and the distribution of the 
clouds may even produce a correlation between R and T (as 
e.g. in the spherical case considered by M01). 

These conclusions should be however tempered by the 
fact that, in Seyfert galaxies, the R-T correlation is still a 
controversial matter on several grounds. Indeed, the measur- 
ment of the reflection amplitude R suffers from large uncer- 
tainties. In particular, it appears that the R values derived 
from spectral fits are very sensitive to the modelling of the 
primary spectrum (see e.g. Petrucci et al. 2001; Perola et 
al. 2002; Malzac & Petrucci 2002) as well as to the eventual 
presence of a secondary components such as a soft excess 
(Matt 2001). Possible effects such as ionization of the reflect- 
ing material, eventual relativistic smearing, Comptonisation 
of the reflected component may all affect the R values de- 
rived from the neutral slab reflection model generally used 
(pexrav). We also note that in NGC 5506, Lamer, Uttley 
& McHardy (2000) find an anti-correlation between R and 

r. 

Also, even if the correlation is real, its interpretation 
could be different. Indeed, following Nandra et al. (2000), 
Malzac & Petrucci (2001, 2002) show that at least part of 
the observed correlation could be understood in term of the 
effects of time-delays between the fluctuations of the pri- 
mary emission and the response of a large-scale distant re- 
flector (e.g a torus). Obviously, because of the time-scales 
involved, this interpretation cannot be valid for black-hole 
binaries. On the other hand, in Seyfert galaxies, it is sup- 
ported by spectral and variability studies indicating that, in 
several sources, an important fraction of the reflected fea- 
tures is produced at large distances from the central engine 
(e.g. Chiang et al. 2000 for NGC 5548; Matt et al. 2001 
and Lamer et al. 2000, for NGC 5506; Done, Madejski & 
Zycki 2000 for IC4329a, Papadakis et al. 2002, for a sample 
of sources) . In a context where the reflection features would 
be produced mainly on distant material, the "cauldron" sce- 
nario considered here would remain a viable possibility. 

Future observational studies of the R-T relation in sam- 
ples of sources as well as in individual objects should clarify 
this point. 
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APPENDIX: FITTING FORMULAE FOR THE 
ESCAPE PROBABILITY 

Let us consider a medium where both scattering and ab- 
sorption of radiation are effective, with Td the total optical 
depth of the medium for the scattering process and r a for ab- 
sorption. Assume one injects instantaneously in the volume 
a number of photons according to a given spatial distribu- 
tion. The photons can disappear from the medium either 
by real escape (i.e. by crossing the volume boundaries) or 
by absorption. Let / be the distribution of the photons at 
their disappearance time, and E the fraction of photons that 
truly escape the medium (i.e. not absorbed). Then, from the 
definition of the escape probability given by equation (^) , it 
is straightforward to show that: 



1_ 
P 



dt 



f(t')dt' 



(19) 



Using a standard Monte-Carlo method, we estimated E 
and the / function for different geometries and source dis- 
tributions. In these calculations the scattering process was 
assumed isotropic and the distribution of scatterers and ab- 
sorbers was homogeneous. We then computed numerically 
P and searched for analytical formulae giving a reasonable 
representation of the numerical results in the range r<j < 10 
and r a < 10. The expressions given below agree within 5% 
with our numerical results. For slab and sphere geometries, 
the optical depths r a and Td are defined along the full height 
of the slab and the radius of the sphere, respectively. 
For a spherical geometry and uniform isotropic injection: 



-r- = 0.75 

p 



1 



0.25r d 



l+0.49r a (l + 0.1r d ) 



1- 



■2.28r a ) ' 16 . (20) 



In the limit r a =0 this formula reduces to that proposed by 
Stern et al. (1995). 

For a spherical geometry and central isotropic injection: 



1 

P 



= {1. +0.45T d [l. + 3.384 10~ 3 T a 168 

exp{3.3(1.09 - 0.10887r a )(10. - r d ) 8 - 
+0.2999r a + (0.346 + O.lTd)^ 681 }]} 
exp (0.944r a + 1.32 10~ 3 r a 2 ) 
1 + 0.56r a ' 



(21) 
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For an infinite slab geometry and homogeneous and isotropic 
source function: 

1 = 1 + 0.392 ^ 183 - 7Ta 



(— 
V1 + 



p Vl + 183.6r, 

Td/3 

(l + 0.08r0- 57 )(l+0.393r a 1 045 ) 



+ 'ILL (22) 



For an infinite slab geometry and photons injected isotropi- 
cally at midplane: 

1 n nno 1 + 123.2r a , 1.211 

P = °- 992 l + 124.9r a eXp[ °- 250r - 

+0.401r d ' 602+(0 - 748Ta -°' 0464rf)0 ' 813 ]. (23) 

These expressions are general and can be applied to any sit- 
uation where both absorption and scattering are important. 
For the situation considered in this work, the scattering r<j 
and absorption r a optical depths depend on the type of ra- 
diation considered (see section ^): 

• For the Comptonised component: Td=Tr and r a =TB 

• For the reflected component: t^—clr^t b , T a =Tr+(l- 
o,r)tb 

• For the soft component: Td =tb, t^—tt 

Note that we neglect the relativistic effects due to the re- 
duction of the Klein-Nishina cross section at high energies. 
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